Effects of thickness in quantum dots at strong magnetic field 
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We study the effects of thickness on the ground states of two-dimensional quantum dots in high 
magnetic fields. To be specific, we assume the thickness to be small so that only the lowest state in 
the corresponding direction is occupied, but which however leads to a modification of the effective 
interaction between the electrons. We find the ground state phase diagram and demonstrate the 
emergence of new phases as the thickness is accounted for. Finally, the wave functional form and 
vortex structure of different phases is analyzed. 



I. INTRODUCTION 

In recent years, there has been an increasing inter- 
est towards two-dimensional quantum dots realized in 
semiconductor heterostructures. Besides applications in 
quantum information and quantum computing, quantum 
dots are interesting in their own right as an example 
of strongly correlated interacting quantum systems, in 
which the strong influence of the magnetic field entails 
prominent diver seness. Moreover, the results extrapo- 
lated from computationally feasible few electron droplets 
are frequently used to understand macroscopic quantum 
phenomena such as the quantum Hall effects. 

In this paper, we report an exact diagonalization study 
of the effects of the layer thickness on the ground states 
of quantum dots in high magnetic fields. Earlier stud- 
ies have pointed outi^ that in a strictly two-dimensional 
parabolic quantum dot in a strong magnetic field and 
with up to seven electrons, the ground states have strong 
correlations that favor either of the two classical configu- 
rations of the types depicted in Fig. [TJ a regular polygon 
with or without an electron at the center. This leads to 
the allowed angular momentum values M = Mmdd + 
fcAM, where k is a non-negative integer, AM = TV or 
TV - 1 (with ft = 1), and M M dd = N(N - l)/2 is 
the minimum angular momentum of the spin-polarized 
quantum dot dictated by the Pauli principle. Alternative 
to the geometric approach is the microscopic composite 
fermion theory (as opposed to the not so succesful mean- 
field composite fermion theory), generality of which also 
extends to larger particle numbers 4 

The vortices in quantum dots with zero thickness have 
been considered in Ref. [4|. Taking into account the thick- 
ness of the quantum dot is, however, occasionally neces- 
sary for understanding of the experimental results 
Since the fractional quantum Hall effect is destroyed in 



FIG. 1: Likely configurations in six electron quantum dots. 



thick systems^, which relates to unbinding transition of 
the vortices at the electron positions^, it is of a more gen- 
eral interest to examine, what type of vortex structures 
are favored at different values of the thickness parameter. 

In Sec. II, we present the microscopic model and a 
method to solve it. Sec. Ill contains the phase diagrams 
and related discussion, followed by an analysis of the mi- 
croscopic states and vortex structures of different phases 
in Sec. IV. Finally, Sec. V summarizes the main results. 



II. MICROSCOPIC MODEL 

Quantum dots formed in the GaAs/Al^Gai-^As het- 
erostructure are modeled for both lateral and vertical 
devices as droplets of electrons in an effectively two- 
dimensional plane confined by a harmonic external po- 
tential. We use an effective-mass Hamiltonian 
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where TV is the number of electrons, m* = 0.067m e is the 
effective mass and A is the vector potential of the homo- 
geneous magnetic field B perpendicular to the quantum 
dot plane. In the calculations, we set the confinement 
strength huoo to 2 meV as its scaling should merely shift 
the ranges of magnetic fields for different phases. For 
convenience, we exp ress leng ths in units o f effective os- 
cillator length I = yjh/m*uj, where uj = \/oOq + (co c /2) 2 
and uo c = eB /m*c is the effective cyclotron frequency. 

To accommodate the thickness of the electron layer, 
we employ an effective interaction potential 
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in which d is comparable to the extent of the wave func- 
tions in the ^-direction and e = 12.7 (CGS) is the dielec- 
tric constant of the GaAs semiconductor medium. Typ- 
ically, an effective interaction is obtained by assuming 
certain form for the potential in the z-direction, solving 
for its ground state, and integrating out the z-direction. 
While the interaction in Eq. ([2|) does not relate to any 
particular wave function, it leads to same qualitative be- 
haviour, especially for smaller as more realistic effec- 
tive interactions and is extensively used in the fields 
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FIG. 2: The ground state angular momentum phase diagram for (a) five and (b) six electrons as a function of the thickness 
parameter d and the magnetic field B. The vertical lines indicate the regular angular momentum difference AM between 
adjacent phases. For six electrons, the angular momentum difference changes depending on the thickness except for the phases 
A and B, as well as for the states with v — 1, |, and |. 



The ground state of Eq. (pQ) is solved by constructing 
the many-body Hamiltonian matrix in the basis of spin- 
polarized lowest Landau level (more accurately the lowest 
Fock-Darwin band since uoq ^ 0) and finding its lowest 
eigenstate by the Lanczos diagonalization. The former 
constitutes a Landau level projection, an approximation 
that is valid at the high magnetic field regime^ In the 
oscillator units, the single-particle wave functions read 

(z\m) = - 7 ±=z m e- z ~ z/2 , m ^ , (3) 

x 1 ' V7rm! v ' 

where z = x + iy. The non-trivial quantities are the in- 
teraction matrix elements (m/, n / |V(ri2)|m, n). Utilizing 
the angular momentum conservation m! + n' = m + n, 
these can be written in terms of 

M l mn = (m + l,n\V(r 12 )\m,n + l) , /,m,n^0, (4) 



III. PHASE DIAGRAMS 

To obtain a phase diagram, we solve the interaction 
energy of the ground state for each angular momentum 
value at each thickness, and determine the angular mo- 
mentum that has the lowest energy for each pair (d, B) on 
a grid. Once we have established the angular momenta of 
neighboring phases, say M\ and M2 for certain d, we ex- 
actly solve the magnetic field B at the phase boundary by 
equating the energies E(d,B,M\) and E(d, B^M^) and 
using a standard routine to solve B. We limit the con- 
sidered angular momenta to M ^ M(y = |) for five and 
six electrons and to M ^ M{y — |) for seven and eight 
electrons, where we utilize the stability of the Laughlin 
states in order to only have stable phases in the phase 
diagrams. The corresponding angular momenta can be 
computed from the trial wave functions^ and are given 

by 



for which an analytic formula is given in the Appendix. 
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FIG. 3: Same as Fig. but for (, 

The single-particle bases are limited according to m ^ 
k(N — 1), for the two fc, which is the basis size needed 
for the y — \ Laughlin state. Returning back to Fig. [TJ 
we remark that according to Eq. ([5]) the states that sup- 
port both classical configurations have the same angular 
momentum as the Laughlin states. 

Fig. [2] presents the phase diagram that is the ground 
state angular momentum as a function of the thickness 
parameter d and the magnetic field B for five and six 
electrons. If the electron number is less than six, the 
structure of the phase diagram remains unchanged as 
the thickness is increased as seen in Fig. [2fa) for N = 5. 
The adjacent phases are then always separated by an- 
gular momentum AM = N. In general, the monotonic 
behaviour of the phase boundaries results from the de- 
crease in the interaction energy as the thickness is in- 
creased. In the case of six electrons in Fig. [2fb), there is 
a transition from AM = N — 1 structure to AM = N 
structure as the thickness is increased. As an exception, 
the MDD state with a hole, M = Mmdd + 6 (phase 
A), has a stable phase all along. Also the phase with 
M = M(y = |) — 6 (phase B), being the analogy of 
the Laughlin's quasielectron wave function, narrowly re- 
tains across the whole thickness range as indicated by 
the break in the vertical line. As a side note, we have es- 
tablished that the overlap between states with the same 
angular momentum (ifj(d = 0)\i/j(d)} decreases monotoni- 
cally but remains significant through the parameter range 
d G [0, 3/] (the overlap is ~ 0.67 for v = | at d = 3/ and 
ranges up to 0.98 for the other states). 

The same basic pattern recurs to the systems with 
seven and eight electrons as seen in Fig. [3j where we 
have confined ourself to the ground states with M ^ 
M(y = |). For eight electrons, there are two additional 
states with M = 46 and M = 52 at small d (shaded 
in the Fig. [31(b)) that appear to be subsequent elements 
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) seven and (b) eight electrons. 

of sequence M — Mmdd + k(N — 2) with k = 3 and 
4. Indeed, we have established that the most probable 
configuration for these states is such that six electrons 
form a ring around a pair of electrons that resides near 
the center of the dot. The latter angular momentum, 
M = 52, is also equal to that of three holes in the maxi- 
mum density droplet M = Mmdd + 3iV (phase C found 
in the diagram after d ~ 2.21) and the Pfaffian stated 
Afpf = Ar ( 2 ^ r ~ 3 ) m j n R e f th e overlap of the corre- 
sponding states with the Pfaffian state was investigated 
and was found to peak high slightly below d = 21 for 
eight electrons. Interestingly, this ground state angular 
momentum is unstable when d G [0.5/, 2.2/] indicating 
that, despite the high overlap, electron pairing is not en- 
ergetically favorable in the lowest Landau level even if 
the thickness is optimized. Finally, we should mention 
the minor phase at d = 3/ and B = 10 T with M = 64 
that can be scarcely seen in Fig. [3](b) and is of little im- 
portance. 

With regard to experiments and an estimate for d, 
we have compared the phase diagram for a two-electron 
quantum dot (figure not shown) to the data presented in 
Ref. where a full 3D treatment of the system was used 
to explain experimental results. The magnetic field val- 
ues at the transition points between the spin-polarized 
ground states appear to be in a good agreement with the 
data when the thickness parameter d ~ 10 nm, which 
corresponds to about d~/atI? = 10T. Therefore, we 
expect at least the effects (besides the obvious shifts in 
B) around d ~ I for six and eight electrons to be relevant 
for an experiment similar to that conducted in Ref. in 
which the electron number, however, was limited to the 
maximum of five. 

At all electron numbers, the range of magnetic field 
for each phase tends to be equal as the thickness is large. 
This behaviour also points towards microscopic explana- 





FIG. 4: The complex phase of the conditional wave functions for the six electron quantum dot, in which the electron co-ordinates 
are fixed according to the most likely configuration indicated by +, and the upmost electron is used as the probe electron. In 
the first row, the thickness parameter d = I and the angular momenta are M = Mmdd + 5/c with k = 2, 3, 4, 5, and 6. In the 
second row, the thickness parameter d — 31 and the angular momenta are M = Mmdd + 6k with k = 1, 2, 3, 4, and 5. 



tion of these phases by vortex formation as the creation 
of an additional vortex is associated with a constant flux 
increase by one flux quantum $0 = hc/ei^ 



IV. VORTEX STRUCTURES 

While the values of angular momenta in the phase 
diagrams are naturally explained by the geometric ar- 
gument, a more accurate description of the microscopic 
states can be obtained by investigating the zeros of the 
wave functions^ Since the many-body wave function is 
up to the exponential factor merely an analytic polyno- 
mial in the complex electron co-ordinates 2^, the zeros are 
clearly visible in the phase structure of the wave function. 
For this purpose, we define the phase of the conditional 
wave function 
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where z[ denote the electron co-ordinates fixed to chosen 
positions. Circumvention of each zero, henceforth called 
a vortex, of p(z) = *&(z, Z2, . . . , zn) accumulates a factor 
2tt. 

The phases cp for the system of six electrons of the 
ground states up to the angular momentum M(y = |) 
at d = I and 3/ are shown in Fig. 0J The vortex on 
top of each fixed co-ordinate follows from the necessary 
presence of the factor Yl( z i ~ z j) * n each solution \£. Of 
special interest are the free vortices in the central area 
of the quantum dot, which typically have an irrelevant 
reflection vortex outside the polygon pattern. In the up- 
per row of the figure, an angular momentum increase by 
AM = 5 always leads to an additional vortex near the 
center, while on the lower row an increase of AM = 6 
is required. Similar results hold for the larger angular 
momenta and electron numbers. 



TABLE I: Overlaps with the trial wave functions for AM = 5 
and AM = 6 vortex states for 6 electron quantum dot. 
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The ground states with M = Mmdd + kN can be 
approximated by the (unnormalized) trial wave func- 
tion $ C M = Y[{zi - Zj)Y[(zi - z C M) k . This is ob- 
tained by eliminating the center-of-mass motion by the 
transformation 15 z\ i— > Z{ — ^cm from the v = 1 integer 
Hall effect state with k holes in the origin. The approx- 
imation scheme is generalized to arbitrary angular mo- 
menta M by taking the highest weight many-body basis 
vector |$) of the ground state and performing the trans- 
formation Zi Zi — ^cm, |$) ► |3>cm)« We denote the 
basis state with quantum numbers mi, . . . , tun by a se- 
quence, in which the elements number mi, . . . , mjsr are 1 
and the rest are zero. The thus obtained trial wave func- 
tions and their overlap with the ground states are pre- 
sented in Table I. At d = Z, the overlaps are rather high 
until M = M(y = |) = 45, which is still better described 
by the corresponding Laughlin state trial wave function 
(overlap ~ 0.97). At d = 3/, the overlaps are higher and 
even M = M(y = |) is slightly better described by our 
trial wave function than as a Laughlin state (~ 0.67). 
The two rightmost diagrams in Fig. [H with d = I (up- 
per) and d = 31 (lower) both having M(y = |) also 
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illustrate this break-down of the incompressible Laugh- 
lin state in thick systems as the vortices are gradually 
less bound to the electron co-ordinates. This effect due 
to the effective interaction is in contrast to the oppo- 
site effect caused by a screened Coulomb interaction^ 
whereby, in the strong screening limit, the zeros are ex- 
actly localized to the electron positions. Both cases are, 
of course, direct consequence of Haldane's pseudopoten- 
tial parameters of the corresponding interactions and the 
fact that the Laughlin state is exactly obtained for an in- 
teraction having only the first pseudopotential coefficient 
nonzero 

The conclusion to be drawn from Table I is that for 
the AM — N the hole is created or extended at the 
origin, while for AM = N — 1 states the essential mech- 
anism of increasing the angular momentum is to create a 
hole next to the origin. Similarly, for the configurations 
with two electrons at the center realized in the eight elec- 
tron dot, the trial wave functions |11000111111)cm and 
|110000111111)cm may be applicable. 



V. SUMMARY 

We have studied the phase diagram of spin-polarized 
quantum dots in high magnetic fields as a function of the 
thickness, and found roughly that for electron number N 
equal to six and seven, the ground state angular momenta 
at different thickness occur at either intervals N or TV — 1. 
For N ^ 5, only the intervals N occur, while for N ^ 



8 other intervals of type N — k can occur. Moreover, 
we have interpreted the angular momenta in terms of 
vortex formation near the center of the quantum dot and 
constructed trial wave functions that have a moderate 
overlap with the exact solutions. 



VI. APPENDIX 

The formula for the interaction matrix elements in the 
effective oscillator units obtained by transforming into 
relative and center-of-mass co-ordinates reads 



M Ln = — [(m + l)\n\m\(n + Z)!]"i 
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where is the step function, \F\ is the Kummer confluent 
hypergeometric function, and I* and a* are the effective 
oscillator length and the effective Bohr radius a* = 47r ^ 2 . 
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